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Abstract. Finite dynamical systems (e.g. Boolean networks and logical models) have been used in 
modeling biological systems to focus attention on the qualitative features of the system, such as the 
wiring diagram. Since the analysis of such systems is hard, it is necessary to focus on subclasses that 
have the properties of being general enough for modeling and simple enough for theoretical analysis. 
In this paper we propose the class of AND-NOT networks for modeling biological systems and show 
that it provides several advantages. Some of the advantages include: Any finite dynamical system 
can be written as an AND-NOT network with similar dynamical properties. There is a one-to-one 
correspondence between AND-NOT networks, their wiring diagrams, and their dynamics. Results about 
AND-NOT networks can be stated at the wiring diagram level without losing any information. Results 
about AND-NOT networks are applicable to any Boolean network. We apply our results to a Boolean 
model of Th-cell differentiation. 



1. Introduction 

Discrete models have a long and successful history in systems biology, beginning with Boolean 
network representations of molecular networks (T| and their later generalization, so-called logical 
models (?). They are qualitative, time-discrete models that are particularly suitable for the analysis 
of steady state behavior of molecular networks. However, as models become larger it is increasingly 
difficult to analyze them. In order to keep the analysis of such networks tractable, many studies have 
focused on specific classes of networks such as: single-switch, unate, nested canalizing, threshold, 
AND, AND-OR, and linear networks (11 [4j [5] [6] [7j [8] [3 [lOj [TT] [T2] [T3] [111 . In order to be useful for 
modeling, a family of networks has to be "sufficiently general" for modeling biological interactions 
and "simple enough" for theoretical analysis. In this paper we propose the family of AND-NOT 
networks as such family. AND-NOT networks are a particular the class of Boolean networks that are 
constructed using only the AND (A) and NOT (-■) operators. 

A biological justification for the use of AND-NOT networks is that there is evidence that for 
genes that are regulated by more than one other gene, the different binding sites exhibit synergistic 
effects between the different regulators |15l|16|[T7l . This fact motivated the study of conjunctive 
Boolean networks, that is, networks whose logical rules are constructed using exclusively the AND 
operator f5j, where explicit formulas for steady states are given; also, upper and lower bounds for the 
number and length of limit cycles are provided. But conjunctive Boolean networks cannot account for 
inhibitory regulation and the resulting negative feedback loops, which are common in gene regulatory 
networks. Allowing the NOT operator, in addition to the AND operator (i.e. using AND-NOT 
networks), can make the family of networks sufficiently general to be useful for modeling 1 1 81 . 

For a formal argument that the family of AND-NOT networks is general enough for modeling, 
we will show that any discrete model (finite dynamical system, to be precise) can be represented 
by an AND-NOT network. More precisely, we present an algorithm that assigns to a given general 
discrete model an AND-NOT network which has the same number of steady states, together with an 
algorithmic correspondence between steady states of the two networks. This is achieved by adding 
nodes to the network as needed. The potential drawback of this algorithm is of course that the 
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network size can potentially get significantly larger, thereby potentially negating any computational 
advantage gained by the specialized logic. However, since molecular networks have typically small 
in-degree, this growth in the number of network nodes to be added is modest in the case of molecular 
network models. We demonstrate this through an analysis of several published models and random 
networks. 

To argue that AND-NOT networks are simple enough for theoretical analysis, we will show how 
using the specialized logic of AND-NOT networks can provide better theoretical results. For exam- 
ple, in (19 1, it was shown that an upper bound for the number of steady states can easily be computed 
for AND-NOT networks (which is not true for arbitrary networks). Also, in 1201 , it was shown that 
the exact number of steady states of AND-NOT networks are encoded in the topological features 
of the wiring diagram, and that, in some cases, the problem of finding the exact number of steady 
states can be transformed to the problem of finding maximal independent sets of the wiring diagram, 
which has been extensively studied 0lll22l[^[24l[25][26l[27j[28][29][30l. In this paper we will show 
how the specialized logic of AND-NOT networks can give us better upper bounds for the number 
of steady states; more precisely, we provide an upper bound for AND-NOT networks that improves 
on previous upper bounds. Furthermore, we show how this upper bound for AND-NOT networks 
can actually be used for general networks. We use our results to analyze a Boolean model of Th-cell 
differentiation. Another theoretical advantage of AND-NOT networks is that they are in a one-to-one 
correspondence with their wiring diagrams. This observation has several implications, one of which 
is the possibility to relate dynamic network properties with features of the wiring diagram (51 1201 . 
Also, from a given signed wiring diagram one can unambiguously construct and AND-NOT network, 
which implies that all algorithms or results can be stated at the "wiring diagram level." 



Definition 2.1. For a signed directed graph G = (Va,Ec), we denote /, = [j : (j,i,s) e Eg], 
= U • C/i ! > + ) e Eg) and IT = [j : (J, i, -) 6 £ G ).That is, /, is the set of all incoming edges for 
node i, and It, resp. It is the subset of positive, resp. negative, edges. All graphs in the rest of the 
paper will be signed directed graphs unless noted otherwise. 

In order to simplify the graphical representation, we denote two negative (positive) edges between 
i and j by a bidirectional negative (positive) edge, •-• (■* — ►). If the edges have different signs we 
denote them by • — ► . 

Definition 2.2. An AND-NOT function is a Boolean function, h : (0, 1)" — > {0, 1), such that h can be 
written in the form 



where P n N = { }. If P = N = j ), then h is the constant function 1. If i e P (i e N, respectively) 
we say that i or x, is a positive (negative) regulator of h or that it is an activator (repressor). An 

AND-NOT network is a Boolean network (BN), / = (/, /„) : {0, 1}" -> (0, 1|", such that f t is 

an AND-NOT function for all i = 1, . . . , n. AND-NOT networks are also called signed conjunctive 
networks. 

Definition 2.3. The wiring diagram of an AND-NOT network is defined by a graph G = (Va,E G ) 
with vertices Vq = {1,. . . ,n} (or {jci, . . ., x„)) and edges E G given as follows: (i, j, +) € E G (((, j, -) e 
Eq, respectively) if x, is a positive (negative, respectively) regulator of /}. Notice that nodes cor- 
responding to constant functions have in-degree zero. Also, the wiring diagram of an AND-NOT 
network contains all the information about the network; that is, we only need to specify the wiring 
diagram in order to define an AND-NOT network. 

Example 2.4. Consider the Boolean network / = (fi,. . .,fe) '■ (0, 1 ) 6 — » {0, 1 j 6 given by 

f(x) = (x 2 A x 4 A — 1JC5 , jcj A x 6 A -1X3 A ->x 5 , 1, x 6 A -1X1 A —1X5, x 6 A -oci, 1). It is easy to see that / is 

an AND-NOT network. Its wiring diagram is shown in Figure[T] 



2. Definitions 
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Figure 1. Wiring diagram of the AND-NOT network in Example 12.41 

As mentioned in the introduction, some other families of networks that have been studied in the 
past are single-switch, linear, AND, AND-OR, unate and nested canalyzing functions (3] [4] [5] [6] 
l7l[8ll9l [T2"l , Each family has its own advantages; however, for the purpose of modeling biological 
systems and for theoretical analysis, it is of interest to have the following properties: First, networks 
generated using these families should be able to admit a sign assignment; that is, it should be possible 
to determine the sign of an interaction. Second, in principle, it should be general enough to model all 
networks; that is, it should be possible to model any type of regulation. Third, for theoretical analysis, 
it would be useful to have a one-to-one correspondence between wiring diagrams and networks. This 
property would allow complete encoding of a network in its wiring diagram. The family of linear 
functions satisfies the third property but not the first two. The family of AND functions satisfies 
the first and third property but not the second. The family of AND-OR functions satisfies the first 
property but not the last two. Single-switch, unate, and nested canalyzing functions satisfy the first 
two properties but not the third. 

On the other hand, AND-NOT networks satisfy all three properties. The first property is satisfied 
because the sign of a regulation is given by the presence or absence of the NOT operator. The third 
property follows from the fact that if the positive and negative edges to i are given by P and N, resp., 
then the function for node i is uniquely given by fi = /\j s pXj A A/ew ~' x i- The second property 
is given by the fact that any finite dynamical system can be expressed as an AND-NOT network. 
More precisely, Theorem l3. 1 I guarantees that steady states are preserved if we rewrite a general finite 
dynamical system as an AND-NOT network. 

3. Results 

In this section we show why AND-NOT networks are a good framework for modeling biological 
systems. 

3.1. AND-NOT networks are general enough for modeling. One issue that can potentially arise 
when only using certain classes of networks is that one can have difficulty in modeling certain pro- 
cesses. For example, the family of AND networks does not allow modeling negative interactions. 
Another example is that the family of linear networks, does not allow modeling signed interactions. 
In order for a family of networks to be useful for modeling, is has to allow modeling any type of 
interaction. 

Here we show that the family of AND-NOT networks is general enough for modeling. More 
precisely, we show that for any finite dynamical system, there exists an AND-NOT network (possibly 
with more nodes) such that they share key dynamical properties. 

Theorem 3.1. Let h = (h u . . . , h„) : S — » S be a finite dynamical system, where S = Xi x ■ • ■ x X„ 
and all Xi 's are finite. Then, there exists an AND-NOT network g : {0, 1 )'" — » {0, 1}"' such that there 
is a bijection between the steady states of h and g. Furthermore, g and the bijection between steady 
states is given algorithmically. We say that g is an AND-NOT representation ofh. 

Proof. A simple proof uses the facts that any finite dynamical system can be written as a Boolean 
network 1311 , and that any Boolean function has a conjunctive normal form. 

In 1311 . the authors proved algorithmically that for any finite dynamical system h, there exists a 
Boolean network / (possibly with more nodes) such that h and / have the same number of steady 
states. Furthermore, the bijection of steady states is also given algorithmically. Therefore, we only 
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need to show that there exists and AND-NOT network g, such that there is a bijection between the 
steady states of / and g. 

We proceed by induction. First, consider the conjunctive normal form oif„: f„ = w\ Aw 2 A- ■ -Aw,, 
where wj is of the form Wj(x) = si X\ V ,? 2 x 2 V • ■ • V s u x„ with j,- s {id, ->} {id =identity function). Notice 
that -iWj is an AND-NOT function. Then, define the BN k = (Jfci... . A+r) : (0, -> {0, lj" +r in 
variables (xi,...,x n ,yi,.. . ,y r ) by £,(x,y) = fi(x) for i = 1, . . . ,n— 1, k„(x,y) = — ■>•] A -iy 2 A ■ ■ ■ A ->y r 
and £,(x,y) = -<Wi(x) for i = n + 1, . . . , n + r. 

We now check that the function (f>{ x) = (x, ->wi (x), . . . , -ivvv(x)) gives a one-to-one correspondence 
between steady states of / and k. Suppose that /(x) = x, then 
k(4>(,x)) = k{f x (x), . . .,jk-i(Jc),Wi(z) A w 2 (x) A ... A w,.(x),wi(x), . . . ,w r (x)) 

= fc(/i(x), ...,/„-i (x),/„(x),wi(x), ...,w r (x)) = (x,-.w,(x), . . . , -.hv(x)) = 0(x); that is, <f>(x) is a 
steady state fc. Now, suppose that fc(x, y) = (x, y) and notice that in this case y; = £;(x, y) = -iw,(x); 
then (x,y) = 0(x). Also, /(x) = (/,(x), • • • ,/„_ 1 (x),/„(x)) = (Mx,y), . . . ,k n „i(x,y), Wi(x) A w 2 (x) A 
••• A w,-(x)) = (x 1 ,...,x„_i,-._vi A -ry 2 A ••• A -._y r ) = , . . . , x n -\ , k n (x, yj) = (xi, . . . , X n -\, X n ) = x. 
That is, x is a steady state of /. Therefore, k = (k[ , . . . , k„ +r ) is a BN where k„, . . . , k„ +r are AND- 
NOT functions and such that there is a one-to-one correspondence between the steady states of / and 
k. By induction, it follows that there is an AND-NOT network g : {0, 1}"' — > {0, 1}"' together with a 
bijection between the steady states of / and g. 

Therefore, there is a bijection between the steady states of h and g. Furthermore, g and the 
bijection are given algorithmically. □ 

The transformation of finite dynamical systems to Boolean networks has been discussed in 1311 . 
So, in the rest of the paper we will focus on Boolean networks and AND-NOT networks. 

Example 3.2. Consider the BN / : (0, l} 5 -» (0, l) 5 given by f { = x 2 V -1X4, f 2 = Xi A x 3 , f 3 = 
(x 2 V -1X4) A x 5 , f 4 = x 3 V x 5 , f 5 = x 3 . The wiring diagram of / is given in Figure [2] (left). In order 
to transform this BN to an AND-NOT network we introduce the variable x 6 with Boolean function 
fb = ~ , x 2 A X4 and fi = -1x3 A -1x5. Variables X6 and X7 will be used in gi and #4. Notice that since 
x 2 V -1X4 appears again in / 3 , we can simply reuse x 6 to keep the number of extra variables as small 
as possible. Then the AND-NOT network is g : (0, l) 7 — > (0, l} 7 given by gi = -iXg, g 2 = X\ A x 3 , 
g3 = -1X6 A X5, g4 = -1x7, g5 = X3, g6 = -ix 2 A X4, g7 = -1x3 A -1x5. The wiring diagram of g is shown 
in Figure|2](right). 



Figure 2. Wiring diagram of the BN network / and the AND-NOT network g in 
Example |3.2| 

An additional step in the transformation that can keep the number of extra variables small is given 
by the following proposition. 

Proposition 3.3. Let / : (0, 1 \" -> {0, 1 )" be a BN and define g : (0, 1 }" -> {0, 1 }" by g = N k ofo N k , 

where Nk(xi, . . . , x„) = (xi, . . . , Xt-i, ->Xt, Xk+i, ■ • • , x„). Then f and g are dynamically equivalent. 

Proof. It is enough to notice that N k is invertible with inverse N k . Then, g r = N k o f o N k ; that is, 
evaluating / is equivalent to evaluating g. □ 

If some functions of a BN are OR-NOT functions, then we can use Proposition [33] to transform 
the BN into a BN in the same number of variables such that the OR-NOT functions become AND- 
NOT functions. Also, Proposition 13.31 can be used to transform constant functions f k = into 
constant functions f k =i (if f k = 0, then the k-th coordinate function of N k o f o N k is the constant 
function 1). 
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Example 3.4. Consider the BN / : (0, l) 3 -> {0, l} 3 given by f\ = x 2 , f 2 = *i V -ue 3 , f 3 = x 2 A x 3 . 
The wiring diagram of / is in Figure [3] (left). Since / 2 is an OR-NOT function, we can transform 
it to a AND-NOT function using Proposition 13.31 Consider g = N 2 ° f ° N 2 , given by g(x) = 

N2(f(Xl,-*X2,X 3 )) = N 2 (^X 2 ,X { V -1*3,-1*2 A X 3 ) = (-1*2, ~i(*l V -1*3), -1*2 A * 3 ) = (-1*2,-1*1 A 

*3, -1x2 A x 3 ), with wiring diagram shown in Figure [3] (right). Then, / is dynamically equivalent to 
an AND-NOT network. Notice that the effect of this transformation on the wiring diagram is simple, 
we simply change the signs of the edges around node 2. 




Figure 3. Wiring diagram of the BN network / and the AND-NOT network g in 
Example [34] 



As mentioned in 1311 , an advantage of transforming finite dynamical systems into Boolean net- 
works is that it can provide insight into the role of feedback loops by disentangling them. In this 
sense, transforming finite dynamical systems into AND-NOT networks can pass all the information 
of the role of feedback loops to the wiring diagram. In this case, the wiring diagram is not only a 
rough representation of the network, but it encodes all the information of the network; in this sense 
the wiring diagram "becomes" the network. This has the potential to reduce the problem of studying 
the structure of the state space graph (which has 2" elements) to studying the structure of the wiring 
diagram of the AND-NOT representation (which has m > n elements). This can help in understand- 
ing the precise role of the network topology in the network dynamics. A similar approach was used 
successfully to study conjunctive and linear networks (5] E) . 

3.2. The variable growth in AND-NOT representation is small. For practical purposes it is im- 
portant to obtain an estimate of how much the AND-NOT representation can increase the number 
of variables. For arbitrary Boolean networks, the number of extra nodes can be exponential in the 
number of nodes. However, Boolean models of biological systems are not arbitrary and are actually 
very sparse with very low in-degree (typically described by a power law distribution 1321 1331 ), We 
will now show that in practice the number of variables introduced by the algorithm can be small. 

Table 1. Number of extra variables introduced by the AND-NOT representation. 
The number of nodes of / and its AND-NOT representation, g, are denoted by 
n, m, respectively. The BNs were taken from |34||35|[36ll37||38| . 



n 


m 


% increase 


12 


13 


8% 


12 


15 


20% 


14 


15 


7% 


20 


24 


20% 


23 


26 


13% 


28 


28 


0% 


40 


43 


7.5% 



In order to study this question, we have applied the procedure to several published models in 
the literature and studied the question using randomly generated Boolean networks. The first study 
shows that the increase in the number of variables for published models is modest (Table [T}. The 
number of variables was increased by 14% on average with a maximum value of 4 extra nodes. In 
order to determine the number of extra nodes introduced by our algorithm for more general BNs, 
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Table 2. Average number of extra variables introduced by the AND-NOT repre- 
sentation for random BNs. 



in-deg< K 


% increase 


K = 1,2 


0% 


K = 3 


5.2% 


K = 4 


10.8% 


K = 5 


16.2% 


K = 6 


20.8% 


K = l 


24.8% 


K = S 


28.6% 


K = 9 


32.3% 


K=10 


36.1% 



we did a statistical analysis. To mimic wiring diagrams coming from biological systems, the edges 
followed a power law distribution and we considered the maximum in-degree less than or equal to K 
for K = 1, . . . , 10 (see Appendix A for details). The results of this second study are shown in Table 
[2] For example, all networks where nodes have in-degree bounded by K = 2 can be transformed 
to AND-NOT networks without increasing the number of nodes. For networks where nodes have 
in-degree bounded by K = 4, our method increases the number of nodes by 10.8% on average (see 
Appendix A for details). It is important to mention that in both tables, the growth in the number of 
extra nodes is far less than exponential. 

3.3. AND-NOT networks can be useful is theoretical analysis. As mentioned in the Introduction, 
the specialized logic of AND-NOT networks can be used to obtain better theoretical results. Such 
results can arise directly (e.g. (51 1201 ) or by applying results about general Boolean networks to the 
family of AND-NOT networks. In this section we show examples of the latter. First, we need the 
following definitions. 

Let C be a feedback loop of a graph G. We say that C is a strong feedback loop if there are no 
edges of the form k — > i,k -» jinG\C such that i, j e C. 

For example, consider the graph G in Figure[4] The feedback loop (3, 4) is not strong because of 
the edges 1 — > 3, 1 -• 4; (5,6) and (1,2,4,3,5) are not strong because of the edges 1 — > 5, 1 -• 4. 
All other feedback loops are strong. 




G 



Figure 4. Graph with only one strong positive feedback loop. 

Our first result in this section is an application of H39I Theorem 3.2] to the family of AND-NOT 
networks (see Appendix B for the proof). 

Theorem 3.5. Let W be the wiring diagram of an AND-NOT network, and suppose J intersects all 
strong positive feedback loops ofW. Then, the number of steady states is at most 2' J '. 

Example 3.6. Consider the AND-NOT network with wiring diagram given in Figure [4] The only 
strong positive feedback loops are (1,2) and (1,3,5). Since J = (1) intersects them, Theorem 13.51 
guarantees that there are at most 2' 7 ' = 2 steady states. 
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Intuitively, Theorem l3.5l is telling us which positive feedback loops contribute to the presence of 
steady states; it says that they have to be strong. We also provide a slight generalization of Theorem 
13.51 We need the following definition. 

A feedback loop C of a graph W is called inconsistent if there is a vertex kc such that there is a 
positive path of the form kc — » f 1 —>■••—>;,—> tc from kc to tc e C and a negative path of the form 
kc ~ * j\ —* • • • — * jr ~* u c> from k c to u c € C such that k c — » fc, fcc — • «c are not edges in C and 
|/,, | = . . . = \Ij r \ = 1. When such vertex kc does not exist, we say that C is consistent. 

For example, consider the graph W in Figure[5] The positive feedback loop (3,4) is inconsistent 
because of the paths 1 — » 3 and 1 — » 2 -• 4. The positive feedback loop {5, 6) is inconsistent because 
of the paths 1 — > 3 — > 5 and 1 -• 6. Also, the positive feedback loop {1,2,4,3,5) is inconsistent 
because of the paths 1 — > 3 — > 5 and 1 — » 2 -• 4. Then, the only consistent feedback loops are (1,2) 
and (1,3, 5). 




w 



Figure 5. Wiring diagram of the AND-NOT network in Example 1 3. 8 1 

We say that a set / C {1, . , . , n] dominates a graph W if J intersects all consistent positive feedback 
loop and for each feedback loop C that is inconsistent and strong, J intersects C or J contains at least 
one kc- For example, the set (1) dominates the graph W in Figure[5] 

With these definitions we have the following theorem that gives an upper bound on the number 
of steady states using topological features of the wiring diagram (see Appendix B for the proof). 

Theorem 3.7. Let W be the wiring diagram of an AND-NOT network, and suppose J dominates W. 
Then the number of steady states is at most 2' J '. 

It is not difficult to see that the bound given by Theorem 13.51 is greater than or equal than the 
bound given by Theorem 1 3. 7 1 The next example shows that the inequality is in some cases strict. 

Example 3.8. Consider the BN / : (0, 1 } 6 -> {0, 1 ) 6 given by 

f = X2 A X$, 

h = x u 

fl = X\ A -1X4, 
/ 4 = -,x 2 A -1X3, 
fs = *3 A ~<X 6 , 

f 6 = -a, a -i^tj, 

Its wiring diagram is shown in Figure|5] It is easy to see that (1,3,5) intersects all strong positive 
feedback loops. Then, Theorem 13.51 gives the upper bound 2 3 = 8. On the other hand, since (1) 
dominates the wiring diagram, Theorem 1 3 . 7 1 gives the upper bound 2. That is, Theorem 13.71 gave 
a better upper bound on the number of steady states. Notice that in this case the actual number of 
steady states is 2, namely, 000101 and 111010. 

One might argue that having better results for AND-NOT networks is not enough to justify their 
use. After all, since we are considering a smaller family of Boolean networks we should of course 
obtain stronger results. However, the combination of Theorem 1 3 . 1 1 and results about AND-NOT net- 
works automatically generates theorems for all Boolean networks. Furthermore, such combination 
can in some cases provide stronger results. This deserves further explanation which is illustrated 
in Figure [6] Consider a theorem about Boolean networks that gives us information about certain 
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dynamical properties, "Thm.". On the other hand, consider a similar theorem about AND-NOT net- 
works, "Thm.*". Then, given a Boolean network /, we have two choices, we can apply Thm. to /; 
or, we can use Theorem 13, II to find the AND-NOT representation of /, then apply Thm.*, and then 
use Theorem [TT] to obtain information about the original Boolean network /. In Section [3~l4l we use 
a published Boolean model to show that the latter can give stronger results. 



Boolean 



Thm. 



Dyn. Prop. 



Thm. 3.1 



Thm. 3.1 



AND-NO T Thrf > Dyn. Prop. 



Figure 6. Extension of theorems about AND-NOT network to all Boolean networks. 



For example, combining Theorem 13.11 and [3~7l we obtain the following theorem. 

Theorem 3.9. Let f be any Boolean network and suppose that J dominates the wiring diagram of 
its AND-NOT representation. Then, f has at most 2' J ' steady states. 

We now show that this theorem can in fact provide a better upper bound for the number of steady 
states. 

3.4. Application to Th-cell differentiation. We apply our results to the BN model proposed in t36l 
for Th-cell differentiation. The model is a BN in 23 variables, / : (0, l) 23 — » {0, l} 23 . Below is the 
list of Boolean functions. The wiring diagram is shown in Figure|7] 



Xi 




GAT A3 


f 


= (a-i VXn) a -ix 22 ; 


X 2 




I FN - B 


h 


= 0; 


X 3 




I FN -BR 


h 


= a: 2 ; 


*4 




I FN -y 


fi 


= (Xl4 V X 16 V X 20 V X 22 ) A 1*19 


A'5 




I FN - jR 


fs 


= X4\ 


«6 




IL - 10 


k 


= Xu 


Xi 




IL - 10R 


h 


= x 6 ; 






IL - 12 


h 


= 0; 


Xg 




IL - 12R 


h 


= x 8 A-.x 2 i; 


Xio 




IL - 18 


fio 


= 0; 


All 




IL - 18R 


fu 


= XioA-ix 2 i; 


Xl2 




IL-4 


fu 


= x\ A -.jc 18 ; 


Xl3 




IL-4R 


fu 


— X[2 A ~~ ] X[-j \ 


A 14 




IRAK 


fu 


= xu; 


A 15 




JAKl 


fis 


= X5 A -1*17; 


A' 16 




NFAT 


fit 


- x 2 i; 


A|7 




SOCS1 


fn 


= x [s Vx 22 - 


A IK 




ST AT I 


fu 


= x 3 Vx ls ; 


Xl9 




STAT3 


fl9 


= xr, 


X20 




ST AT 4 


fx 


= Xg A — iJC] ; 


A'21 




STAT6 


/21 


= A43; 


a; 2 




T -bet 


fn 


= (xu V X 22 ) A -iXu 


A'23 




TCR 


fa 


= 0. 



Using our algorithms we obtain the AND-NOT network, g : {0, l} 26 — » {0, l) 26 , shown in Figure 
[8] It turns out that the set (1, 22} dominates the wiring diagram of g (see Appendix C for details). 
Then, by Theorem 13.91 the number of steady states of / is at most 2 2 = 4. On the other hand, 
all previous results about steady states (e.g. 181 1191 ) give 8 as the upper bound. That is, using the 
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Figure 7. Wiring diagram of the Th-cell differentiation model. 



AND-NOT representation can provide a better upper bound, even for general Boolean networks. The 
actual number of steady states of the model is 3 (see 1361 for details). 
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25 26* 
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Figure 8. Wiring diagram of the AND-NOT representation of the Th-cell differ- 
entiation model. 



4. Discussion 

The results presented in this paper, together with other results in the literature, support that the 
family of AND-NOT networks are general enough for modeling and simple enough for theoretical 
analysis. Given any finite dynamical system, it is possible to create an AND-NOT network such that 
they have similar dynamical properties. This has two implications: First, this means that using AND- 
NOT networks in modeling does not pose any technical restriction on the type of interactions one can 
model. Second, every result about AND-NOT networks can be applied to general Boolean networks, 
which can give better results (e.g. Theorem |3.9l >. One potential drawback for this framework is that 
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the AND-NOT representation can have more nodes. However, for networks that arise from modeling 
biological systems, this increase in the number of nodes is modest (Section [3~2ji. 

Other advantages of using AND-NOT networks are the following: First, all information about 
the network is actually contained in the network's wiring diagram. Specifically, there is a one-to-one 
correspondence between AND-NOT networks and graphs, so that the network can be reconstructed 
unambiguously from the wiring diagram. In 1401 the authors followed a similar approach to suc- 
cessfully study cascading effects. Second, due to this correspondence, we can state all results about 
AND-NOT networks using wiring diagrams only. This means that questions about AND-NOT net- 
works can be reformulated as questions about graphs; then, one can use tools from graph theory and 
combinatorics to study them (e.g. antichains, posets, inclusion-exclusion principle, independent sets 
|5, 20]). This deserves further investigation. 

Finally, we point out that AND-NOT networks are special cases of so-called nested ccmalyzing 
Boolean networks. These were first introduced in 110111 II as good candidates for models with "bio- 
logically meaningful" regulatory rules, and have since been studied extensively. In (131 this concept 
was generalized to multi-state models, and it was shown there that the large majority of regulatory 
rules that appear in published models of biological networks are of this form. It was shown further- 
more that nested canalyzing networks have dynamic properties one would expect to find in biological 
networks, such as short limit cycles and a small number of attractors. Thus, the results in the present 
paper imply that in order to study the steady state behavior of general network models, one can focus 
on the very restrictive class of nested canalyzing networks 1141 . instantiated as AND-NOT networks 
and make use of their very special properties. 



Appendix A 

We describe here the details of the study to determine how many nodes are added by the construc- 
tion of the AND-NOT representation. To mimic wiring diagrams coming from biological systems, 
the edges followed a power law distribution. More precisely, given K fixed and a parameter y, the 
probability for a node to have k < K nodes is p k = lc y (up to a normalization factor). For example, 
if K = 4, the probabilities of having 1, 2, 3 and 4 nodes are p\ = cl~ r = c, p 2 = c2~ r , pi = c3~ y and 
P4 = c4~ r , respectively, where c = i- 7+2 - 7 ' +3 - r+ 4-7 so mat Pi + Pi + Pi + P4 = 1. Also, to mimic bi- 
ological regulation, we restricted our analysis to Boolean functions that admitted a sign assignment 
for the edges. These Boolean functions are called unate, biologically meaningful and regulatory 
functions (7]|9][8). 

Denote with e k the average number of extra nodes introduced by a Boolean function in k variables. 
Then, a BN that follows the distribution mentioned above will have, on average, p\e\ + p 2 e 2 + • • ■ + 
Pk£k extra nodes. Now, we need to estimate e k . 

Consider a Boolean function, h, that depends on k variables. For k = 1 there are 2 functions, 
h = x\ and h = -ixj and we do not need to introduce any new nodes; then e\ = 0. For k = 2 there 
are 8 functions and they are of the form h = SjX t A s 2 x 2 or h = 5|Xi V s 2 x 2 , where .s,x, = x, or 
SjXj = -<Xj. For functions of the form h = s t Xi A s 2 x 2 we do not introduce any new nodes, and for 
functions of the form h = sixi V s 2 x 2 we can use Proposition [33] to transform h to an AND-NOT 
function, so we do not introduce new nodes either. Then e 2 = 0. For k = 3, there are 72 functions. 
An exhaustive-search analysis shows that of those 72 Boolean functions, 16 introduce nodes, 48 
introduce 1 node, and 8 introduce 3 nodes; then the average number of extra nodes in this case is 
g3 _ 16*0+4^.1+8.3 _ j p Qr ^ _ ^ there are 1824 Boolean functions. An exhaustive-search analysis 
shows that of those 1824 functions, 32 introduce nodes, 320 introduce 1 node, 480 introduce 2 
nodes, 960 introduce 3 nodes and 32 introduce 4 nodes; thus the average number of extra nodes in 
this case is e 4 = 32«o+320,i+480.^2 + 960»3+32*4 _ 235. For k = 5, there are 220608 functions and an 
exhaustive-search analysis shows that e 5 = 4.03. For k = 6 there are approximately 5 x 10 8 functions 
and an exhaustive-search analysis would be unfeasible. However, we have the following result. 

Theorem A.l. The average number of extra nodes for a unate function ofk variables is at most 
C(k, L&/2J); that is, e k < C{k, |_fc/2J). Where C is the binomial coefficient and [ J is the floor function. 
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Proof. Without loss of generality we assume the CNF of the Boolean function / has no negative 
signs. Let / = vfi A . . . A w r be the CNF, where w-, has the form Wi = X\ V . . . V x s . For each i, define 
St = {/ : Xi appears in w,j. 

Now, if there are i, j such that 5; C S j, then we can simplify w, A wj to w, (e.g. (xi V x%) A {x\ V 
x 2 V X3) = x\ V X2). That is, we can simplify the CNF so that S-,%S j for all i t j. 

Thus, 5 1 , . . . , S r is a family of subsets of j 1 , . . . , k j such that no one is contained in the other. 
Sperner's theorem (4TJ states that r < C(k, [k/2]). This implies that for any unate function in k vari- 
ables, we need at most < C(k, [k/2j) extra nodes to obtain the AND-NOT representation. Therefore, 
e k < C(k, L*/2J). 

□ 

It is important to mention that the exhaustive-search analysis done for k = 3, 4, 5 suggests that e k 
is actually much smaller than C(k, lk/2j). In fact, we did a statistical analysis for k = 6, . . . , 10 using 
a total of 5000000 Boolean functions chosen at random (1000000 for each k). The analysis shows 
the following approximations: e 6 w 5.32, e 7 « 7.04, eg « 9.32, e 9 « 12.24, e l0 w 15.96. 

Table|2]shows a summary of our analysis for y = 2.5. For example, if K = 4, then the fractions 
of functions with 1, 2, 3 and 4 variables are on average pi = .786, p 2 = .139, p 3 = .0504 and 
p 4 = .0246, respectively. Then, the average number of extra nodes is: 

100(p,e! + p 2 e 2 + p 3 e 3 + p 4 e 4 ) = 100(.786 * + .139 * + .0504 * 1 + .0246 * 2.35) ~ 10.8%. 

Appendix B 

Here we prove Theorem 1 3 . 5 1 and [3771 As mentioned in Section [3~3l Theorem |3.5| is an application 
of 1391 Theorem 3.2] to the family of AND-NOT networks. First we need the following definition. 

Let / : (0, 1)" (0, 1}" be a Boolean network and consider x e (0, lj". Then, W(x) = (V, E) is 
the graph with vertices V = {1, . . . ,n] and the following edges: 
0', +) e E if Xj = and / ; (x) < fi(x + ej), or if Xj = 1 and fi(x - e,) < / ; (x); 
(j, i, -) € E if xj = and f,(x) > fj(x + ej), or if Xj = 1 and f,{x - ej) > fj(x); 

where ej is the vector given by (e^), = Sjj (6 is the Kronecker delta). Notice that if (j, i, +) or (J, i, -) 
is an edge in W(x), then changing the j-th coordinate of j produces a change in fj. Notice that for 
AND-NOT networks we have that W(x) C W for all x; in fact, this is true for more general networks. 

Theorem B.1. I39I Let f be a Boolean network and suppose a and b are steady states of f. Then, 
there there exists x such that W(x) has a positive feedback loop with vertices in the set {i : Oj ^ bj\. 

We now prove Theorem 1 3. 5 1 

Proof. Let <f> : (0, 1)" — > {0, l} 1 '" defined by <p(x) = xj. We will show that if a j= b are steady states 
of g, then <p(a) t <f>(b). Consider a + b steady states of g; then, by Theorem B.I., there exists x such 
that W(x) has a positive feedback loop, C, with vertices in the set {i : a, j= b,}. 

We claim that C is a strong positive feedback loop of W. By contradiction, suppose there is 
k e {1, . . , , n] and i, j e C such that k — » i and k -» j are edges in W(x) but not in C. Then, W(x) has 
edges of the form (Zi, ±) and (l 2 , j, ±) where li,l 2 ± k. On the other hand, since (k, i, +), (k, j, -) e 
C c W(x) c W, we have that / ; = x k A . . . and f = -<x k A . . .. We have two cases x k = or x k = 1. 
In the case x k = we obtain that /; = for all values of jci, . . . , in, Jft+i, . . . , x„. In particular, W(x) 
cannot have an edge of the form (/, i, ±) with I t k; this is a contradiction. In the case x k = 1 we 
obtain that /, = for all values of Xi , . . , , Xk-i , JCt+i , . . . , x n . In particular, W(x) cannot have an edge 
of the form (/, /', ±) with / # k\ this is a contradiction as well. Therefore, C is strong. 

Since C is a strong positive feedback loop in W, C must intersect J. Since C has all its vertices 
in the set {( : a, t £>,}, J intersects the set {( : a, t Z>,). Therefore (f>(a) = aj + bj = <f>{b). 
It follows that the restriction of (f> to the set of steady states is an injective function. Therefore, 
|{x:/(x) = x)|<|{0,l} w | = 2'-'l. □ 

It is important to mention that Theorem 1 3. 5 1 was also proven in 1191 using different techniques. 
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We now prove Theorem 1 3. 7 1 

Proof. Let / : {0, lj" -» {0, 1|" be an AND-NOT network with wiring diagram W. Let C be a 
positive feedback loop that is strong and inconsistent. Then, there is a vertex k c such that there is a 
positive path of the form kc — » i\ —»•••—» i r — > ?c from fee to fc e C and a negative path of the form 
kc — > y'l — > ■ ■ ■ — > J, ~~* «c> from fe c to u c e C such that k c ^ t c , k c -• Uc are not edges in C and 
|/„ | = . . . = \Ij r |=1. Let G be the graph obtained by adding to W all edges of the form kc — » 
and £ c -• Uc where C does not intersect J. Denote by g : {0, 1}" — » {0, 1}" the AND-NOT network 
associated to Z. We claim that the steady states of / and g are the same. We prove this by induction 
on the number of extra edges. 

Suppose that W and Z only differ in the edge k — * t, then, by definition we must also have a path 
k — * i\ —»...—» — > /. Suppose that g(x) = x, we need to show that fj(x) = Xj for all j. Since 
W and Z only differ in the edge k — > ? we have fj = gj for j + t, g, = f, A jc; and /, = x, r A . . .. 
Then, fj(x) = g ; (x) = x ; - for j t. It remains to show that f,(x) = x t . Consider first the case x, = 0, 
then, g,(x) = and xt = for some i e /, + . If i 4- k, we have that the edge i — > ? is in W and 
/< = Xi r A Xj A . . .; then, / f (x) = x ir A A . . . = = x t . If i = k, then X* = which implies that = 
(because of the edge k — » ii); similarly, we obtain that x ir = 0. Then, /,(x) = A . . . = = x t . That 
is, f,(x) = x,. Now consider the case x, = 1. Since 1 = x, = g t (x) = f t (x) Ax*, we have /,(x) = 1 = x,. 
A similar argument shows that if f(x) = x, then g(x) = x. The proof for when W and Z only differ in 
the edge k -» t is analogous. By induction we obtain that / and the AND-NOT network obtained by 
a completion of W have the same steady states. 

Now, we claim that J intersects all strong positive feedback loops of Z. Let C be a strong positive 
feedback loop of Z. Then we have two cases: C is in W or it is not. Consider the case C C W. 
Then, C is a strong positive feedback loop in W. If C is consistent in W, then it intersects J. If C 
is inconsistent (and strong) in W, then it also intersects J. Now consider the case C $£ IV. Then, at 
least one edge of C is of the form k c — > /c or A c -• k c for some C strong and inconsistent that does 
not intersect 7. Then, kc E J and 7 intersects C. In any case we obtain that J intersects all strong 
positive feedback loops of Z. 

Then, the number of steady states of g, and hence /, is at most 2 |y| . □ 

Appendix C 

We first analyze the original BN using previous results. In 1361 , the authors showed that the 
positive feedback loops of the BN / : {0, 1 1 23 -> / : {0, 1 ) 23 are: 
{4,5, 15, 18, 12, 13,21, 11, 14) 
{4,5,15,18,12,13,21,9,20) 
{4,5,15,18,12,13,21,1,6,7,19) 
{4,5, 15, 18, 12, 13,21, 1,20) 
{4,5, 15, 18, 12, 13,21, 1,22) 
{4,5,15,18,17,13,21,11,14) 
{4,5,15,18,17,13,21,9,20) 
{4,5,15,18,17,13,21,1,6,7,19) 
{4,5,15,18,17,13,21,1,20) 
{4,5, 15, 18, 17, 13,21, 1,22) 
{4,5,15,18,22) 

{4,5, 15, 18,22, 17, 13,21, 11, 14) 
{4,5, 15, 18,22, 17, 13,21,9,20) 
{4, 5, 15, 18, 22, 17, 13, 21, 1, 6, 7, 19) 
{4,5, 15, 18,22, 17, 13,21, 1,20) 
{4,5, 15, 18,22, 1, 12, 13,21, 11, 14) 
{4,5, 15, 18,22, 1, 12, 13,21,9,20} 
{4,5, 15, 18,22, 1,6,7, 19) 
{4,5, 15, 18,22, 1,20) 
{12, 13,21, 1) 
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(13,21, 1,22, 17) 

(22) 
(22, 1) 
U) 

We will use the following two theorems (proven in f8l U9l , respectively) that give upper bounds 
on the number of steady states. 

Theorem 4.1. Let W be the wiring diagram of a BN network and suppose J is a set of vertices that 
intersects all positive feedback loops in W. Then, the number of steady states is at most 2 |J| . 

Theorem 4.2. Let W be the wiring diagram of a BN network and suppose J is a set of vertices that 
intersects all functional positive feedback loops in W. Then, the number of steady states is at most 
2 M 

It is easy to see that all positive feedback loops intersect the set (1,4,22). Therefore, Theorem 
14. II gives the upper bound 2 3 = 8. Also, it is possible to show that the functional positive feedback 
loops are (4, 5, 15, 18, 12, 13, 21, 11, 14), (22), (22, 1) and (1) (e.g. using the GLNsim software (42) ). 
Therefore, Theorem |4.2| gives the upper bound 8 as well. 

We now analyze the AND-NOT network using our results. The positive feedback loops of the 
AND-NOT network in Figure[8]are the following (new nodes are in bold). 

{24,4,5, 15, 18, 12, 13,21, 11, 14) 

(24,4,5, 15, 18, 12, 13,21,9, 20) 

(24, 4, 5, 15, 18, 12, 13, 21, 26, 1, 6, 7, 19) 

(24, 4, 5, 15, 18, 12, 13, 21, 26, 1, 20) 

{24, 4, 5, 15, 18, 12, 13, 21, 26, 1, 22) 

{24,4,5,15,18,17,13,21,11,14) 

(24,4,5, 15, 18, 17, 13,21,9, 20) 

(24, 4, 5, 15, 18, 17, 13, 21, 26, 1, 6, 7, 19) 

(24, 4, 5, 15, 18, 17, 13, 21, 26, 1, 20) 

{24, 4, 5, 15, 18, 17, 13, 21, 26, 1, 22) 

{24, 4,5,15,18,25, 22) 

{24, 4, 5, 15, 18, 25, 22, 17, 13, 21, 1 1, 14) 

{24, 4, 5, 15, 18, 25, 22, 17, 13, 21, 9, 20) 

{24, 4, 5, 15, 18, 25, 22, 17, 13, 21, 26, 1, 6, 7, 19) 

(24, 4, 5, 15, 18, 25, 22, 17, 13, 21, 26, 1, 20) 

(24, 4, 5, 15, 18, 25, 22, 1, 12, 13, 21, 1 1, 14) 

(24, 4, 5, 15, 18, 25, 22, 1, 12, 13, 21, 9, 20) 

{24, 4, 5, 15, 18, 25, 22, 1, 6, 7, 19) 

{24,4,5, 15, 18,25,22, 1,20) 

(12, 13,21, 1) 

(13,21,26, 1,22, 17) 

(22, 25) 

(22,1) 

{1,26) 

Those feedback loops that contain 4 and 13 are inconsistent because of the paths 1 — * 12 — » 13, 
1 — > 6 — > 7 — » 19 — • 4; they are also strong. All other positive feedback loops are consistent 
and intersect (1,22). That is, (1,22) intersects all consistent positive feedback loops, and for each 
positive feedback loop C that is inconsistent and strong, / contains k c = I. Hence, (1, 22) dominates 
the wiring diagram of g. Therefore, Theorem |3. 91 gives the better upper bound 2 2 = 4 on the number 
of steady states of /. 
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